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We consider the tails of probability density functions (PDF) for different characteristics of ve¬ 
locity that satisfies Burgers equation driven by a large-scale force. The saddle-point approximation 
is employed in the path integral so that the calculation of the PDF tails boils down to finding 
the special field-force configuration (instanton) that realizes the extremum of probability. We cal¬ 
culate high moments of the velocity gradient d x u and find out that they correspond to the PDF 
with \n[P(d x u)\ oc — {—d x u/ Re) 3 / 2 where Re is the Reynolds number. That stretched exponen¬ 
tial form is valid for negative d x u with the modulus much larger than its root-mean-square (rms) 
value. The respective tail of PDF for negative velocity differences w is steeper than Gaussian, 
ln'P(ry) ~ -(w/iirms) 3 , as well as single-point velocity PDF In"P(u) ~ — (|u|/u rms ) 3 . For high 
velocity derivatives u ^ = d x u, the general formula is found: ln'Pdw^l) oc —(|u ( - fc ^|/Re i: ) 3/, ^ +1 b 


1. INTRODUCTION 


The forced Burgers equation 

dtu + ud x u — vd x u = (f> ( 1 - 1 ) 

describes the evolution of weak one-dimensional acoustic perturbations in the reference frame moving with the sound 
velocity [Q. It is natural to assume the external force </> to be ^-correlated in time in that frame: 

( <j)(ti,xi)(t){t 2 ,X 2 )) = 5(ti - t 2 )x(x i - x 2 ). ( 1 . 2 ) 

Then the statistics of (j> is Gaussian and is entirely characterized by 0 )- We are interested in turbulence excited by 
a large-scale pumping with the typical correlation length L so that x does not essentially change at x < L and goes 
to zero where x > L. Besides L , the correlation function x is characterized by the parameter to = [— (l/2)x"(0)] 1 / 3 
having the dimensionality of frequency. Then e.g. x(0) ~ L 2 lo 3 . Developed turbulence corresponds to the large value 
of Reynolds number Re = L 2 lo/v. 

The physical picture of Burgers turbulence is quite clear: arbitrary localized perturbation evolves into shock wave 
with the viscous width of the front, that gives k~ 2 for the energy spectrum at Re kL 1 [|l]|^]. The presence 
of shocks leads to a strong intermittency, PDF of velocity gradients is substantially non-Gaussian || and there is 
an extreme anomalous scaling for the moments of velocity differences w = u(p) — u(—p ): ( w n ) oc ( p/L ) for n > 1 
||]. Simplicity of the equation and transparency of underlying physics make it reasonable to hope that consistent 
formalism for the description of intermittency could be developed starting from Burgers equation 10 as well as from 
the problem of passive scalar advection p fll|. That would be an appropriate way to mark centenary of Burgers’ 
birth. 

The most striking manifestation of intermittency is the strong deviation of the PDF tails from Gaussian or the 
behavior of high moments. The PDF’s V(d x u) and V(w) are not symmetric, the asymmetry is due to the simple 
fact that positive gradients are smeared while the steepening of negative gradients could be stopped only by viscosity. 
It has been realized recently that those tails can be found by considering the saddle-point configurations (we call 
them instantons) in the path integral determining the PDF ]l2| |. The instanton formalism has been applied to the 
Burgers equation by Gurarie and Migdal [^3| who found the right tail InV(w) ~ — [w/(ujp)] 3 determined by practically 
inviscid behavior of smooth ramps between the shocks. That cubic right tail has been first predicted by Polyakov 
from the conjecture on the operator product expansion it corresponds to the same right tail for the gradients 
lnV(d x u) ~ — (d x u/to) 3 which also has been derived from the mapping closure by Gotoh and Kraichnan 0- Here 
we find the viscous instantons that give the left tails of PDF’s and single-point velocity PDF. 

Even though some calculations are lengthy, the simple picture appears as a result. Since white forcing pumps 
velocity by the law w 2 oc (f> 2 t while the typical time of growth is restricted by the breaking time t ~ L/w then the 
Gaussianity of the forcing In?'’(<(>) oc —<(> 2 /x(0) leads to \nV(w) ~ —[|iu| /{Lu>)] 3 . At a shock, w 2 ~ —ud x u so that 
lnV(d x u) ~ -[-9i;ii/(wRe)] 3 ^. These simple estimates are confirmed below by consistent calculations. 
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2. SADDLE-POINT APPROXIMATION 


A general instanton formalism for the description of the high-order correlation functions and of the tails of probability 
density functions has been developed in [|l2| . Here, we will be interested in the high-order moments of the velocity 
gradient d x u , difference w = u{x = p) — u{x = —p) and the velocity itself: 


{[d x u] n ) = Z 0 1 J VuVpexp{iX + nh\[d x u(Q, 0)]} , 

(w n ) = Zq 1 j VuVpexp {il + nln[w(0, p) — u(0, — p)]} , 


( 2 . 1 ) 

( 2 . 2 ) 


= z; 


'Du'Dpexp {il + nlnit(0,0)} . 


(2.3) 


Here, Zq is the normalization constant and the effective action X = f dtC(u,p). The form of the Lagrangian C for 
(1.1 |l.2| ) is determined by the dynamical equations |[fo] [iT| : 

(2.4) 


C = J dx (pd t u + pud x u - vpdlu) + ^ J dxidx 2 piXi 2 P 2 ■ 


The expressions (2.1,2.2) with (2.4) are formally exact relations. 


A. Saddle-Point Equations 


The main idea implemented here is that the h igh-ord er correlation functions (for n 1) are determined by the 
saddle-point configurations of the path integrals ([0,2.2). The corresponding saddle-point equations are 


We should substitute here 


d t u + ud x u — vd x u = —i J dx 1 x( x ~ x ')p( x ')i 

(2.5) 

d t p + ud x p + vdIp = 5(t)X(x). 

(2.6) 

A => [in/a(Q)\5'(x) a(0) = d x u(0,0) for gradients, 

(2.7) 

A => —(in/w)[6(x — p) — S(x + p)} for differences. 

(2.8) 

A => — [m/u(0,0)]<5(cc) for velocity. 

(2.9) 


A solution of the equations ( 2.5 2d]) determines the saddle-point value of the arguments of the exponent in i |2. l[ffo] i 
and consequently the corresponding correlation functions in the steepest descend approximation. 

The eq uati on (2.5) can be solved with arbitrary initial condition at remote negative time. The condition for the 
equation (2d>) should be formulated for the final time [^2|. It is p = 0 as it follows from the variation of the action (|2.4| ) 
over u at the final time, since the corresponding contribution to the action o rigin ating from f dtdxpdtu is f dxpSu. 
That means that p = 0 at t > 0. Thus we can solve the system of equations (2.5) and 


d t p + ud x p + v dip = 0, 


( 2 . 10 ) 


at negative time t only. The role of the last term in (2.6) is then reduced to the initial (or more precisely final) 
condition p{x) = — A(x) imposed on p at t = — 0. The initial condition imposed on it at a remote negative time is 
forgotten after a finite time and does not influence the answer. From the other hand the final condition on the field p 
is also forgotten at a finite negative time if to move backwards in time. Thus our solution is characterized by a long 
stationary stage where both fields u and p are zer o. W e will refer to that stage as vacuum since it is realized at the 
absence of the source A in the right-hand side of (243). In the following we will treat the vacuum value u = 0 as the 
initial condition for the field it. As to the field p it grows from a small perturbation at moving in the positive direction 
in time due to the negative sign at the viscous term in ( 2.10| ). We can say that our instanton solution describes the 
instability inheren t to the vacuum. 

The equations (2.5 2.1C ) are invariant under x —> —x, u —> —u, p —> —p. Thus for gradients and differences, the 
antisymmetry of the boundary term A leads to an antisymmetric solution u(x) = —u{—x), p(x) = —p(—x). For 
velocity, the solution will be asymmetric. 
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In the saddle-point approximation 


([ d x u) n ) ~ exp(*X extr )[a(0)] n , ( w n ) ~ exp(zJ e xtr)[w(0)]” , (u n ) ~ exp(zT extr )[it(0, 0)] T 


( 2 . 11 ) 


where we should substitute the corresponding values found as a solution of (2^,16) and the saddle-point value of 
the effecti ve actio n. Of course, besides the saddle-point contributions, there exist fluctuation corrections to the path 
integrals (EM ( 3 ]). The corrections at the vacuum stage are compensated by the normalization constant Zq since 
the vacuum just corresponds to n = 0 that is to A = 0. The contribution of the fluctuations around the instanton 
(associated with nonzero saddle-point values of the fields u and p ) is finite and small in comparison with the saddle- 
point one if n 3> 1. 


B. Conservation Laws 


The equations (2.5, 2.10 ) are Lagrangian equations and the Lagrangian ( p2. l| ) does not explicitly depend on time. 
Then the “energy” E is conserved which can be found by the standard canonical transformation: 


iE = / dtpdtu — C. 


( 2 . 12 ) 


The explicit expression for E can be obtained from (2.4): 

E = i J dx (jpud x u + vd x pd x u) - (1/2) J dx\ dx -2 piXizP'i ■ 


(2.13) 


Since the Lagrangian (2A) does not explicitly depend on coordinates then the “momentum” J is conserved as well: 


iJ = 


dxpd x u. 


(2.14) 


Because of the conservation laws we should treat solutions of (2.5 2.6) with zero values of both E and J since they 
are zero in the va cuum . Using E = 0 we conclude from (2.4 2.13| ,2.12) that the saddle-point value of the effective 
action figuring in ( 2.11 ) is 


T e xtr — / dt dx pd t u. 


(2.15) 


Any conservation law is an important heuristic tool to be used at the beginning to understand general properties of 
the solutions. Let us consider the final instant t = —0. Here p(x) = —A(x). Now, substituting that into the energy, we 
can analyze the balance of different terms. For the gradient, E = —na(0)+uj 3 n 2 /a 2 (0)—ni'd x u(0, 0)/a(0) = 0. Striking 
difference between the cases a(0) > 0 and a(0) < 0 is now clearly seen from the energy conservation. Indeed, for the 
case of positive gradient, the viscous contribution to the energy is unessential and two first terms can compensate 
each other. Then one gets a(0) = wn 1 / 3 , which corresponds to the right PDF tail In V{d x u) ~ —(d x u/u>) 3 given by 
an inviscid instanton (see the next subsection). The instanton that gives a(0) < 0 cannot exist without viscosity as 
it is seen from energy conservation. The same consideration can be developed for velocity differences with the same 
result that positive w comes from an inviscid instanton while negative from a viscous one. For the velocity, we use 
both conservation laws. Requiring J to be zero we obtain d x u(0 , 0) = 0. That leads to the following energy at t = —0: 
E = x(0)n 2 /w 2 (0) — z/<9 2 it(0,0). Without viscous term, energy conservation cannot be satisfied. Let us emphasize 
that the answer we shall obtain for the tails of velocity PDF (as well as that for velocity differences) does not contain 
viscosity while it’s consistent derivation requires the account of the viscous terms in the equations. 


C. Right Tail 

Let us first describe the instantons producing the right tails of the PDFs for gradients and differences f| |l3| , |l4[ . 
The tails correspond to the positive slope of the velocity u near the origin. Then one can expect that the field u is 
smooth enough and the viscosity plays a minor role in its dynamics. Due to the initial condition p = —A at t = 0 the 
field p is localized near the origin. At moving backward in time the viscosity will spread the field p. Nevertheless, 
if u > 0 at x > 0 and u < 0 at x < 0 then the sweeping will tend to “compress” the field p. In this situation one 
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can expect that the width of p remains much smaller than L. Then, it is possible to formulate the closed system of 
equations for the quantities a(t) = d x u(t , 0) and 


=(t) = -i/<ix*p(t,*), 


(2.16) 


regarding that the width of p is so s mall t hat the velocity u can be replaced by it’s linear expansion term ax. The 
equation for c c an be obtained from (2.1C) if to multiply it by x and integrate over x. The equation for a can be 
obtained from (225) where for narrow p we can put f dx'x(x — x')p(t,x r ) —> — id x x{x)c(t ) and for small x we can 
substitute d x x{x) « —2 uj 3 x. We thus find 


d t c = 2 ac, d t a = — a 2 + 2 w 3 c. 


(2.17) 


The energy conserved by (2.17) can be derived from ( 2.13| ) at the same approximation: E = —ca 2 + ui 3 c 2 . The 
condition E = 0 means that the instanton is a separatrix solution of ( 2.17 ). 

For gradients p(—0 } x) = — [in / a(0)\8' (x) i.e. a(0)c(0) = n. From the energy conservation w 3 c 2 = ca 2 , one gets 
a(0) = w 3 c 2 (0)/n = a m 1 / 3 . The same solution describes also the instanton for differences since in this case p(— 0, x) 
is determined by (|2.8|) what leads to w = 2a(0)p and then to the same initial condition a(0)c(0) = n. Note that the 


life time of the instanton is T ~ n 1 . 


The main term in ( 2.11 ) is [a(0)] n , that immediately leads to (( d x u) n ) ~ [a(0)] n ~ w n n ra / 3 which gives the right 
cubic tail of the PDF In V(d x u) ~ — (9 x u/u>) 3 . Analogousl y on e can obtain In V(w) ~ — [w/(pu>)] 3 found previously 
JuJEl. Then one can check that the extremum value (p.!5|) is negligible in comparison with nln[a(0)]: 


2extr = i J dtcdta ~ a(0)c(0) = n. 


One can incorporate viscosity and show that it has negligible influence on the answer and that the width of p is 
much less than L through the relevant time of evolution giving the main contribution into the action ]l3[ . The right 
tails of V(d x u) and V{w) are thus universal i.e. independent of the large-scale properties of the pumping. 


D. Left Tail 


The main subject of this paper is the analysis of the instantons that give the tails of Viu) and the left tails of 
V(d x u) and V(w) corresponding to n egative values of a(0), w. Although the held p is narrow at t = 0, we cannot use 
the simple system of equations (2.17) to describe this instanton as it has been seen from the conservation laws. In 
terms of the function u, the reason is in the shock forming near the origin which cannot be described in terms of the 
inviscid equations. In addition, sweeping by negative velocity slope provides for stretching (rather than comp ress ion) 
of the held p at moving backwards in time. Thus we should use the system of full differential equations ( |2.5[ . [f~fi| ,i 
to describe the instantons which can be called viscous instantons. Before presenting analytical consideration, let us 
qualitatively describe the instantons and make simple estimates of their parameters. 

We shall see that, apart from a narrow front near x = 0, the velocity held has L as the only characteristic scale of 
change. The life time |Tj of the instanton is then determined by the moment when the width of p reaches L due to 
sweeping by the velocity uq : |Tj ~ L/uq. Such a velocity uq itself has been created during the life time |Tj by the 
forcing so that rto ~ |c| max TLo; 3 . To estimate the maximal value of |c(t)|, let us consider the backward evolution from 
t = 0. We first notice that the width of p (which was zero at t = 0) is getting larger than the width of the velocity 
front ~ uo/a already after the short time ~ a -1 . After that time, the values of c and a are of order of their values at 
t = 0. Then, one may consider that p(t,x) propagates (backwards in time) in the almost homogeneous velocity field 
Mo so that dtc = — if dxxup x ~ uo / dxp. The integral / dxp can be estimated by it’s value at t = 0 which is ti/uq. 
Therefore, we get c max — nT so that |Tj ~ n” 1 / 3 ^ -1 and zto — Lwn 1 / 3 . 

Since at the viscosity-balanced shock, the velocity u o a nd th e gradient a are related by Mq — va then a(0) ~ 
wRen 2 / 3 . The main contribution to the saddle-point value Q is again related to the term [d x u(Q, 0)]" and we find 
((d x u) n ) ~ [a(0)] n ~ (wRe)™n 2n / 3 , which corresponds to the following left tail of PDF at d x u wRe 


V{d x u) oc exp[— C{— 9 x u/wRe) 3 / 2 ] . (2-18) 

The velocity difference w{p) is uq for L^> p^> x o where the width of the shock xq — uo/a(0) ~ n _1 / 3 Re _1 . We thus 
have ( w n ) ~ ( Lu>) n n n / 3 which corresponds to the cubic left tail 

V{w) oc exp{— B[w /(Leo)] 3 } . 
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The product Luj plays the role of zt rms . The instanton describing the velocity statistics is asymmetric since the velocity 
maximum should reach x = 0 at t = 0 as it has been seen from the momentum conservation. Nevertheless, all above 
estimates concerning w are valid also for u and we obtain 

V(u) oc exp{— D[u/(L lo)] 3 } . 

The numerical factors C, B and D are of order unity, they are determined by the evolution at t ~ T i.e. by the 
behavior of pumping correlation function xi x ) at x ~ L. The left tails are thus nonuniversal; in addition to L and ui, 
they are determined by the large-scale properties of the pumping. 


3. ANALYTIC DESCRIPTION OF THE VISCOUS INSTANTON 


The formation of the shock ne ar the or igin occurs at time £ 
inviscid version of the equations ([T5 2.1(j|) 


< 


\T\ and at small x. At moderate x we can use the 


dtu + ud x u=—i j dx'x{% ~ x')p(x'), 
d t p + ud x p = 0 . 

That system obeys a rescaling symmetry which admits the following form of solutions 

u = Nf 1 (t/T,x/L ), p=N 2 f 2 (t/T,x/L), T=-N~ 1 uj- 1 , 


(3.1) 

(3.2) 

(3.3) 


where fi, f 2 satisfy corresponding equations. We will demonstrate that a unique family with N rel ated to n is 
realized for our viscous instanton. To prove that, we should find a solution of the viscous equations (2.5, 2.1fj| ) at 
i<L, |f| -C |T| and to match the solution with (3.3) at intermediate scales and time. 

We shall construct this short-scale solution at this section. Here, we want to note only that the shock formed at 
small time near the origin c an b e in the main approximation described by the stationary equation ud x u — vd 2 x u = 0 
which can be derived from (2.5) where the time derivative and the pumping term can be neglected in comparison 
with large space derivatives. Then we find the well known expression 


2v 

u = — tanh(ay:co), 
Xo 


(3.4) 


describing the field u at t = 0 and sma ll enough x, the width of the shock xq should be mu ch le ss than L. At scales 
x xq we return to the behavior (3.3) what gives the estimate LujN for the amplitude of (3.4). Then we find from 

(PD 


Xo 


LRe _1 £V _1 , a(0) ~ w7V 2 Re, w ~ LujN, 


(3.5) 


where w = w(0, p ) — u(0, —p) for xo p <C L. 

Since the instanton describing the velocity statistics is asymmetric the shock is not formed near the origin. It can 
be described by the expression u = ( 2 v / xq ) tanh[(x — p)/xo] + const obtaining from ( jOl ) by shifting. Thus the general 
conclusions made above are valid for the instanton but the velocity statistics needs a separate consideration. The 
next subsections are devoted to the instantons describing the statistics of the derivatives and of the differences and 
the last subsection presents the instanton for the velocity statistics. 


A. Cole-Hopf Representation 


It is known that the equation (1.1) can be rewritten in the linear form by using Cole-Hopf substitution [|TJ. To 
proceed further we shall introduce the substitution for our instanton. Namely we reformulate the problem for the 
new variables and P: 


= 2 vd x p=P<Z. 

2v 


(3.6) 


The action (2.4) is rewritten as 
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1= / dt dxP(d t ^ — vd 2 ^>) + - / dtdx\dx2Xi2PiP2- 


( 3 . 7 ) 


The energy ( 2.1 3| ) is now 


E = i / dx vd x Pd x ty + 


iv 


dx FP^ 


= —iy dxPdt^f — j dxFPty, where 

d * F = ~~^2 Jdx’xlx-x’)p{x'). 


(3.8) 

(3.9) 


We will fix F by the condition F (x = 0) = 0. Note that for the antisymmetric solution (describing the gradients and 
the differences) 


F{x = +oo) — F{x = —oo) = — 


2v 2 


dxx{x) / dx'p(x') = 0. 


The saddle-point equations for the functions T and P read 


dt’fy - vd 2 ^ + vF f = 0, 

d t P + vd 2 x P - vFP - 2v\\x)S{t)^- 1 = 0. 


(3.10) 

(3.11) 


Here, we have introduced the term from (2.6), producing the final condition for P, into the corresponding equation. 
One concludes from (3.10|,3.11) that 


<9 t (PT) + vd x (d x P^ - Pd x T) = 2vS(t)X'(x), 


(3.12) 


and the quantity f cfeP’F is conserved. Actually, the integral is equal to zero since P(x) = 0 at t = +0. T hat means 
that we can treat functions p{x) tending t o zer o at x —> ±oo since / dx PT oc J dx d x p in accordance with (3.6). Note 
that the saddle-point value of the action (2.1 5|) can be rewritten also in terms of the fields and P: 


lextr = / dtdxPdt 'F. 


(3.13) 


As we already pointed out, the spatial width of p(t, x) is small at small |i|. It enab les one to go further in formulating 
the equations for the fields P and T. That allows us to expand xi x ~ x ’) i n (3-9) in the series over x'. The leading 
term gives 


F =^M0)-X(x)]. 

Here c is defined by ( 2.1C|) , c can be also rewritten in terms of the fields and P 

c = — / dxx 2 Pty. 

dv J 

For x L ( 3.14 ) is reduced to F(x) = [w 3 /(2 v 2 )\cx 2 . Then, the system ( |3.1C , 3.11 ) can be rewritten as follows 

d t ^> — vd 2 \ P + — cx 2 ^ = 0, 

LO 3 

d t P + vd 2 x P - t^cx 2 P - 2v6(t)\ , (x)<I>- 1 = 0, 


(3.14) 

(3.15) 

3 

(3.16) 

(3.17) 


and the energy (|3.8| ) is 


E=-i dx Pd^ - to 3 c 2 . 


(3.18) 


6 























B. Generating Function 


The temporal evolution of p for time where it’s width is less than L can be effectively described in terms of the 
generating function for the moments 


Y(rj) — —i J dx exp(— rjx 2 /2)d x p . 


(3.19) 


Note that Y(0) =0 and 9^(0) = c where c is determined by the expression (2.16). We shall express the behavior 
of the moments of p vi a some de rivatives of u taken at the initial time t = 0. Since the field p is narrow enough at 
small time we can use ( 1.16 |3.17 ) at x <C L to calculate the generating function ( 3.19| ). It is useful to pass to the new 
variables 


1 

v = —= exp 

VA 


(-£*■)*• 


1 

P = —= exp 

VA 




with the parameters that satisfy ordinary differential equations 

9 t A = ipA, d t <p + (p 2 — 2 uj 3 c = 0. 


(3.20) 


(3.21) 


Since the quantities A and tp are introduced by the differential equations where initial conditions are not fixed, then 
we can impose two arbitrary conditions without changing any observable quantity. Note that the equation (3.21) for 
pi describes the evolution of the linear profile u = px of the velocity which the pumping c would produce at i«l 
in accordance with the inviscid equations. 

Let us introduce the dimensionless variables 


y = x/A, ds = — 2vA 2 dt. 


Then the system (3.16,1.17) leads to the following equations for the functions W and P 

ds'L + = 0, 

d s P - ^d 2 P =-2u5(s)d y X i?- 1 . 


(3.22) 

(3.23) 

(3.24) 


Now, we see the essence of the transfor mation (3.20,1.21): the equations (3.23,1.24) are separated and do not contain 
the pumping explicitly. The expression S) is now rewritten as 


Y = -— J dyP(s,y)ty(s,y)enp(-rjy 2 /2), 

where fj = A 2 rj. The energy ( |3.18 ) in terms of the new variables takes the form 

2 iv 


E = 


A 2 


dyPds'ty + tip / dy y P dy^ + Lo a c — ap . 


The formal solution of the equations ( 3.23| .3.24) can be presented as 

&(s,y) = exp[—s9^/2]T(0, y), 

P{s,y) = —2z/exp[s9y/2]4 , ~ 1 (0, y)d v \. 


(3.25) 


(3.26) 


(3.27) 

(3.28) 


The expression ( 3.25|) becomes Y = if dy ^(s, y) exp(— rjy 2 /2) exp[s9^/2]4' 1 (0, y)d y X. Integrating over dy by parts 
and substituting ( (1.27 ) we obtain 


Y = i I dyty 1 (0, y)d y Xexp •( —\[y + sSy] 2 ! ^(0, y). 




Here we used the following relation 


exp(s9^/2) exp {-py 2 /2) exp(-s9^/2) = exp -My + sd v \ 


(3.29) 


(3.30) 
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At t = 0 near the origin there is the shock (3.4). Taking into account also the linear profile we find at i«i 


't'(C), ir) oc cosh 


x 0 


exp 


1 

4z/ 


(fo + ViW 


(3.31) 


where fo is the value of if at t = 0. We have introduced also the parameter ifi whic h ac counts for the de viatio n of 
the actual solution from tanh(x/iro) at t = 0, f\X is the linear pr ofile to be added to (3.4). Starting from ( 3.31 ) one 
can find the explicit expression for the generating function Y via ( 3.29 ). First we use the Hubbard-Stratonovich trick 
to convert the second-order operator at the exponent into the first-order one: 


exp exp + itlv + sd y } 


d£ 


1 


1 


exp ' ) exp (~ n? s ) exp ^ y ) ex p (*£ sd v ) • 


(3.32) 


We used the relation exp(t\y + t 2 d y ) = exp(iii 2 /2) exp(tiy) exp (t 2 d y ). Substituting ( 3.31 . 3.32 ) into ( |3.29 ) we obtain 
after integrating over £ 

f} (l + (p 2 s) 2 y 2 + s 2 /y$ ~ 

2 1 + fjs(l + (p 2 s) 


Y = i I dydyX 


\/l + 7?S(1 + f>2 S) 


exp 


cosh 


y 


1 + 77 s (1 + (p 2 s) yo 


cosh 


-1 V 


Vo 


(3.33) 


where y 0 = x 0 /A 0 and (p 2 = Ao(^ 0 + y>i)/(2i/) 


C. A solution 


To analyze the t-dependence of (3.33) we should know a solution of the system ( [3,21 ) to restore s(t) from ( |3.22 ). 
For this, it is convenient to use another set of dimensionless variables 


2v _ 2u> _ 2v _ 2u> _ . r /, r 

i P=-^ ( P = — l P’ a= 1 ^ 2 a= — a ' A = L(Re) 7 A. 


A 2 


A 2 


A 2 


Then the equations (3.21) are rewritten as 


d s A = —if A, d s f + (p 2 + -A 6 c = 0, uc = A 2 c. 

Let us write also the equation for the time t following from ( 3.22| ) 

d s t=-— A 2 . 

2u 

One derives from ( |3.20 3.23 , 3.3]]) the following relations 

1 
2 


d s ln'b(s 1 y = 0 ) = ^(a - <p) => a 0 = -1/yl - <p\ 


where = A 2 f\/2v. After expanding (|.3.33|) we find 


iv, 


c = ——(1 + f 2 s) / dydyX 


y z ( 1 + (p 2 s) + 2s— tanh — 


y 


(3.34) 


(3.35) 


(3.36) 


(3.37) 


(3.38) 


y 0 y o_ 

We see that c is the function of the second order over s. Substituting now the formal solution ( 3.27 , 3.28 ) into ( 3.2C| ) 
we find after some calculations using the same tricks as above 


— = A 4 c 2 — Ac(p 2 /A 2 + 4:(if/A 2 )d s c — (2/A 2 )d 2 c. 
u> 


(3.39) 


One can directly check that the energy ( [3.39| ) is c onserved as a consequence of (3.35, [3.3(| ). No w, to establish the 
f-dependence of c we should solve the system ( 3.35 ) and then reconstruct the function s[t) from ( 3.36 ). 

We will see that the last term in the equation ( [3.35 ) for Cp can be neglected in the wide interval of changing s. 
Inside that interval (to be specified below), the equation ( 3.35 ) has a solution 


ip = z 


A = Aqz 




A 2 

ut = --2-Slfo, 
2 z 


(3.40) 


implying that <fo < 0. We see from ( 3.40j ) that A grows with decreasing z and consequently the behavior ( 3.40| ) is 
destroyed for small enough z. 
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D. Gradients 


For the gradients, A is determined by (^dj). Thus one should substitute A = in(2vdo) 1 5'{y) into ( 3.33| ): 

1 


Y _ n (1 + sip 2 )fj f 1 I2s + s 2 (l + sip 2 )fj\ ~ . 

2vdo [1 + s(l + s<p 2 )y} 3 / 2 1 Vo [1 + s(l + s^ 2 )??] ^ ' P 


s 2 fj 


2yl 1 + S7j(l + (p 2 s)_ 


Expanding (3.41) one finds 


c= — {l + sip 2 ) 
2a 0 


2s _ . 

— + (1 + sip 2 ) 
IV 0 


(3.41) 


(3.42) 


It is time to recall that the energy should be equal to zero at o ur so lution. At t = 0 the value of A is very small and 
it is possible to neglect the first term in the right-hand side of ( 3.39 ). Then, equating the residue to zero (and taking 
s = 0) we find <p 2 ~ <fo that is we can believe ipi = 0. That means that with our accuracy it is possible to ne glect the 
linear profile in the the velocity at t = 0 in comparison with the shock contribution. Neglecting (p\ in ( j3.37 ) we find 
also dg = — 1/j/q. Substituting the relations into (3.42) one obtains 


c = mpoz 


-s + —<A) z 
2ao 


(3.43) 


where 2 is determined by ( |3.40 ). We see that c tends to zero where 2 —> 0. It is related to growing A with decreasing 
1 2 1 since the first term in the expression ( 3.39| ) fo r the energy E cannot tend to infinity. 

For 2 <C 1 we can disregard the second term in (3.43) and write c = — nsipoZ. At 2 -C 1, it follows from the definition 
of 2 that s ss —ifQ 1 and we conclude that c « nz. The behavior (3.40) is observed if 1 2 1 3> z\ where z\ ~ ti 1/,3 Aq since 

< 21 , both A and (p do not 


at 2 ~ — 21 the term proportional to A 6 in (3.35) is getting substantial. In the region \z 
vary essentially and can be estimated as follows 


A 


' 1 , 3 Ap, v~n-'' 3 Ap. 


Using fl3.35|]3.40| ) one obtains for 21 <C \z\ <C 1 


c = 2 nt. 


(3.44) 


(3.45) 


This r esult is correct, in particular, for the time where the width of p is much lar ger t han the width of the shock . 
Then ( |3.45 ) has to be compared with c[t) = N 2 fsjt/T) what is the consequence of (3.3) and of the definition ( [2.16 ). 
We thus conclude N = n 1 ^ 3 . Then we find from (3.5) 


*( 0 ) 


r 2/3 Re, T = -wn~ 1/3 . 


(3.46) 


It will be convenient for us to no rmal ize A such that A = L at 2 = 0, then from ( [3.34 ) it follows that A = (Re) 1 / 2 
there. Comparing this value with ([3.44) we conclude 


An ~ Re -1 / 2 n -1 / 3 


Zl 


1//3 Re 


-1 


(3.47) 


Then we find from ( 3.34 . 3.40| ) do 


AgRen 2 / 3 


(p q. It can be rewritten as yoPo ~ 1- Then the second contribution 


to c in (3.43|) is ~ nz 2 and is consequently much smaller than the first contribution at small 2 what .justifi es it s 


neglecting. Now we can find the width of the field p at t = T. We see that at 2 —> 0 all coefficients at 77 in ( 3.41 ) 
tends to zero excepting for the coefficient in the argument of the exponent which for small 2 is exp[— s 2 ^/(2 j/q)]. Just 
the term determines the width of the field p at small 2 , the width is sA/yo ~ A since at small 2 s ~ — (p^ 1 and 
y 0 (po ~ —1. In accordance with (3.34,3.40) the width of the field p increases with |t|: A ~ n 1 / 3 Lut at |t| <C |T| and 
reaches the value A ~ L at t ~ T. 


Now we can estimate corrections associated with neglecting the first term in the expression for the energy ( 3.39| ) 
at t = 0. Using ( 3.42 ) we obtain from (|3.39|) at t = 0 


E 


—A 4 - 

4fig 0 


2 n 


A 2 ho \Vo 


~2 $1 + Vl 


where ipi = <p >2 — <Po determines the linear profile of the velocit y. Th e term with tpf can be neglected here and we find 
ipi ~ nAgipQ. Let us estimate the correction to the expression (3.43) for c following at nonzero ip± from (|3.42|). First, 
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due to the presence of <pi, 1 /j/q differs from —do, see ( 3.37 ). The ratio ipi/ao ~ A® is negligible, see ( 3.47 ). Second, 
the term 1 + sip 2 is now ipo(z + ipiipg 2 ). The question is in the value of the correction at z = z 1 , it is estimated 


as ipi <fi n 2 /zi 


R e~ 2 n~ 2 '\ 


see (3.47), the correction is negligible. Thus we proved the correctness of the expression 
( 3.43 ). Note that the anomalously small value of the linear profile p\ shows that the formation of the shock is finished 
at t = 0 since ip\ directly determines the time derivative of the width of the shock: dt In Xq = ipi ■ This is not surprising 
taking into account extremal properties of the instanton. 


E. Velocity differences 

Let us apply the same scheme for the differences. We introduce the dimensionless quantities 


A 0 w 


P = 


P 


2v ' Ao 

so that we find from (2.8) A = in[6(y + p) — 5(y — p)]/( 2vw). Substituting this expression into ( |3.33 ) we obtain 

fj[(l + sip 2 )p + s/y 0 } 21 


(3.48) 


v _ n [(1 + stp 2 )p+ s/yo] 

+ ^ [1 + s(l + sp 2 )fi} 3 / 2 XP 


2[1 + s(l + sp 2 )fj] 


(3.49) 


At calculating ( 3.49| ) we substituted cosh by exp because of the inequality p yo that means that the separation p 
between the p oints is much larger than the (viscous) width of the front. 

Expanding (3.49) one finds 


c = — (1 + sp 2 ) 


(1 + sp 2 )p + 


Vo] 


(3.50) 


Again, using the relation E = 0 with the energy ([3.39 ) we find that <fi = 0. Then disregarding <pi also in ( 3.3l| ) we 
derive the relation w = —2/yo. Finally, 


c = ni ^° (—2 ws + pipoz). 
w 


(3.51) 


Now, we should s olve t he equations ( 3.35 ) with c given by ( 3.51 ). Again, we have the be havio r ( |3.40 ) stopped by 
the term with c in (|3.35| ). To estimate corresponding z\ we should use the asymptotic of ( 3.51 ) at small z. As we 
shall see, the leading term is the linear one so that c ~ nz. Then we return to the same situation as for gradients. 
That gives 


Lwn 1 / 3 , w ~ Re 1 / 2 A 0 n 1//3 , 


(3.52) 


where we used the condition p xo to be satisfied. To justify the above assertions we should check that at z ~ —Z\ 
the second term in (3.51) is small in comparison with the first one. Utilizing p ~ Re 1//2 p/(LAo) and (3.47) we obtain 
the estimates for the ratios of the noted terms 


P 

—z 1 
w 


LA 2 


luTzi 


P «-V3 


< 1 , 


(3.53) 


Thus the term is small since p/L <C 1 is implied. All the other estimates are the same as for the case of gradients. 


F. Higher derivatives 


The above formalism allows us also t o find t he PDF tails for the odd derivatives u ^ of the ve locity . The aver¬ 
age value ([u/ fc ^] ra ) can be written like (2.1,2.2) what leads to t he saddle-point e xpressio n of the (2.11) type. The 


corresponding instanton is determined by the equations (|2.5|.2.6) where instead of (2 1 7,2 1 8) one should take A to be 


A = 


i( fe )(Q,0) 


S (k \x). 


(3.54) 


All above expressions containing A will be correct for odd k. Substituting (3.54) into (3.38) one finds 
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k + 1 


-nsipoz 


(3.55) 


Performing the analysis similar to that of subsection |3 D| we obtain that N ~ [(k + 1 )n] 1 / 3 and therefore the life 
time T of the instanton is T ~ w -1 [n(fc + l)] -1 / 3 . It becomes smaller for derivatives of higher order and therefore 
our approximation works better for larger k. Then from (|3.4|,3.5) one finds the characteristic value u^ fc ^(0,0) ~ 
N k+1 L 1 ~ k ujRe k leading to 


<J[u (fc) ] n ^) ~ wRe fc L 1 - fc n( fc+1) / 3 . 


The result (3.56) can be rewritten in terms of PDF: 


(V fc) i) 


oc exp 


-c k 


iW\L k 


_1 v 3/(fc+l) 


wRe 


(3.56) 


(3.57) 


Note that the non-Gaussianity increases with increasing k. On the other hand, the higher k the more distant is the 
validity region of ( [1.57 ): u^ 3> uSs ~ L 1 ~ k u)Re k . 


G. Extremum action 

Let us now calculate the extremum value of the action ( 3.13| ). For s < l^ol” 1 one can rewrite it using (3.21,3.21) 
jT extr = — i j dsdy Pd s ^> + J dsc{2p 2 — cA 6 ) — i J dsdyyPdy^p. (3.58) 

Using now the expression ([3.26) for the energy E we can rewrite (3.58) as 


tie. 


= -\f ds (-^ E ~r 2A6+d , 


(3.59) 


Let us recall that the energy E in (3.59) is 0. Now we can estimate all the terms in ([3.59| ) utilizing the asymptotics 
( 3.40 ) valid at 2 > z\ where z\ ~ n 1 / 3 Aq. Then f ds A 2 ~ n -1 / 3 . Substituting ( [j.5C| ) we find f dsc 2 A 6 ~ n. The 
last contribution in (3.59) can be rewritten as 


J dsd = J dta « —21n[4'(0,0)/4'(T,0)], 


the last relation is the consequence of (34^, 3.10). We see that only one contribution to i- T extr is proportional to n, it 
is ~ n. Thus with our accuracy we can neglect the contribution. 

It is natural that p-dependence of V(w) cannot be found in a saddle-point approximation; as a predexponent, it 
can be obtained only at the next step by calculating the contribution of fluctuations around the instanton solution. 
This is consistent with the known fact that the scaling exponent is n -independent for n > 1: (w n (p)) oc p. All said 
above is applicable also for high derivatives. 


H. Statistics of velocities. 

Let us now consider the statistics of velocity u. Equations for fields P and 4/ are the same as for odd derivatives: 

dt® - vd 2 A> + vFV = 0, (3.60) 

d t P + vdlP - vFP - 2v\'{x)8(t)'H- 1 = 0. (3.61) 

while A = —inS(x)/u( 0,0) is not an odd function of x so that F is different now. Since f pdx is nonzero then F is a 
linear (rather than quadratic) function of x for narrow p: 

F(x) = b = —L (dxp(x) (3.62) 

v- 2 J 
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This substantially simplifies the calculations. Substitution which separates variables 


^ = A exp P = A exp (~ipx)P(t,y), 


(3.63) 


includes y = x + p, A and <p which satisfy the following equations 

d t ip + ^—^-b = 0, dtp — = 0, d t A + v<p 2 A = 0. 


(3.64) 


The equation for A is separated and since A does not enter any observable we can forget about it. The equations on 
and P are 


d t ^! - vd 2 ^ = 0, d t P - vd 2 P = -^^-S'(y - p 0 )S(t)^ 1 

wo 


Formal solutions of ( 3.65 ) are 

2 ivn 


P = 


wo 


exp(-vtd 2 )S'{y - p 0 )§ 1 (0,y), § = exp(i/f<9)))'F(0, y) 


Substituting ( 3.66 ) into b = i f dyy^P/Av we obtain 


(3.65) 


(3.66) 


b = —- 


2w 0 


dyS’(y - p 0 )^> 1 {0,y)exp(-vtd%)yV(t,y) 


= y-2vt" v 


d v ^{Q,y) 

*(0 , y ) 


The conservation of the momentum (2.14) gives 

2 ivn 


iJ= dxPdx'S = / dyPdy'i’ = 


U 0 


dyS'(y - p 0 ) 


d v ^(0,y) _ _in 
*{0,y) w 0 


dx5'(x) [u(0, x) + tp] = 0, 


(3.67) 


(3.68) 


what leads to d x u(0, 0) = 0. That natura lly m eans that the maximum of u reaches the point x = 0 at t = 0. We see 
that J is the same integral which enters (3.67). Therefore, b is a constant (given by it’s value at t = —0): 


b= 2PP 0 j dy ^ V ~ = 


2u(0,0) 


(3.69) 


Similarly to the consideration at Sec. p D| , we can now estimate the value of u(0, 0). Velocity stretches the field p 
so that the width of p reaches L at T ~ L/u(0,0) while the velocity itself is produced by the pumping during the 
same time: it(0,0) ~ y(0 )bT = x(0)nT/2u(0,0) ~ nx(0)L/ u(0, 0 ). That gives u(0, 0) ~ Leon 1 / 3 . As well as in the 
preceding subsection, we can neglect the term with X extr in (2.11) and conclude that 

V(u) ~ exp[— D(u/Llu) 3 ] . 

Here D is some numerical factor of order unity which is determined by the evolution at t ~ T i.e. by the large-scale 
behavior of the pumping. 


4. CONCLUSION 

At smooth almost inviscid ramps, velocity differences and gradients are positive and linearly related w(p) = u(p) — 
u(—p) « 2 d x up so that the right tails of PDF have the same form 

V{d x u) ~ exp[-(9 x u/o;) 3 ], V(w) ~ exp[-(u>/2 pu) 3 ] . 

Those tails are universal i.e. they are determined by a single characteristics of the pumping correlation function y(r), 
namely, by it’s second derivative at zero ui = [— (l/2)\"(0)] 1 / 3 . Contrary, the left tails contain nonuniversal constant 
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which depends on a large-scale behavior of the pumping. The left tails come from shock fronts where w 2 — —vd x u so 
that cubic tail for velocity differences corresponds to semi-cubic tail for gradients: 

V{8 x u) ~ exp[— C(—d x u/ojRe) 3 ^ 2 , V(w) ~ exp[— B(w/Luj) 3 ] . 

The formula for V(u>) is valid for w u lms where zt rms ~ Lu>. It is natural to assume that the probability is small 
for both u(p) and u{—p) being large simultaneously. Therefore, V(w) should coincide there with a single-point V(u). 
Indeed, we saw that at u it rms 


V(u) ~ exp[— D(u/Luj) 3 } . 


Here, B,C and D are (related) nonuniversal constants dependent upon the behavior of y(r) at r ~ L. Note that 
the cubic tail for the single-point PDF has been obtained for decaying turbulence by a similar method employing 
the saddle-point approximation in the path integral with time as large parameter 119 . Also, our formula for the 
derivatives V^u^l) oc exp[— Ck(\u^ |i fe_1 /wRe fe ) 3 ^ fe+1 ^] coincides with that of |l9j] for the particular case of white 
(in space) initial conditions. That, probably, means that white-in-time forcing corresponds to white-in-space initial 
conditions. One should be cautious, however, comparing the results for forced and decaying turbulence. 

Another important restriction is our assumption on delta-correlated pumping. If the pumping has a finite correlation 
time r then our results, strictly speaking, are valid for u, w -C L/t. 
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APPENDIX A: 

Consistent account of fluctuations and calculation of predexponent in PDF will be the subject of future work. Yet 
one may worry if it is possible that fluctuations destroy instanton i.e. give some divergent contribution into the action. 
That may happen due to some hidden symmetries and particular fluctuations that are deformations along symmetry 
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group. A special feature of our problem is associated with d-correlated character of the pumping which leads to a 
quasi-symmetry of the gauge-invariance type. Namely, under transformations 


the action X = J dt C determined by Lagrangian 
is equal to zero for functions p satisfying 


p(t, x) —> p(t, x — r), u(t,x)—* u(t,x — r) + dtr, r = r(t ), (Al) 

is transformed as X —> X + f dtdxpdfr. This additional term 


dl / dxp = 0. 


(A2) 


For antisymmetric objects [gradients u' and differences w with A being S'(x) and 8{x + p) — 8(x — p) respectively], 
X(x) = —X(—x). Then the integral J dtdxXu entering ( |2.1 2.2) is invariant under the transformation (|Al]) if we 
require r(t = 0) = 0 which will be implied below. 

We thus have to do with the quasi-gauge transformation which left the effective action invariant on a wide subclass 
of functions. That means that the integration over the quasi-gauge degree of freedom should be performed exactly. 
The conventional way to integrate over the gauge degree of freedom is to pass to the action with a fixed ga uge, 
excluding the volume of the corresponding gauge group. To perform the program for the transformation (Al) we 
include into the integrand in (|2.1[|2.2|) the additional factor 


Vr S[u(t , r) — d t r\J , J = det 


S[u(t, r) - d t r\ 


Sr 


(A3) 


where the integration is performed over functions r(f). The factor (^d|) is equal to unity by its construction and 
therefore d oes not influence Z. Let us now change the order of integration to J Vr J VuVp , and perform the trans¬ 
formation & Then the effective action X is shifted and 6[u(t, r) — d t r} —> S[u(t, 0)]. The Jacobian J introduced by 
(A3) is reduced to 


J -> det[-5* + d x ii(t, x)| x=0 ]. 
Then the integration over r gives the functional ^-function 


Vrexp ( i / dtdxd^rp 


dx d^p{t, x) 


(A4) 


(AS) 


imposing the constraint on the field p. 

Now, we should find the Jacobian (A4) which depends on the regular ization of the transformation ( |Al|) . R eally, 
this regulariza tion is fixed by the regularization of the integral (l j2.l||2.2| ). The point is that at deriving ( |2.4| ) from 
the equation (0) the retarded regularization was implied (otherwise, some additional u-dependent term appears 
originating from the corresponding Jacobian [p~8| ). The retarded regularization means that at discretizing time we 
should substitute 


o , o . u n u n—l / 

O t U + U0 x U -b Un-lU^n 


(A6) 


where e is the step in time and u' = d x u. The discretized version of the transformation & is 

u n (x) -> u n (x - r n ) + fn, 


(A7) 


where f = dtr and r n is the expression to be determined. Substituting (|A7|) into the right-hand side of (A6) we obtain 


u n u n — l / If/ \ f Ml / I • / I 

---h u n -iu n _ 1 — [u n (x - r n ) - u n -i(x - r n -i)\ + +r n -iu n _ 1 +r 

+ ^[u n (x - r„) - u n -i(x - r„_i)] - ^(r n - r n - 1 )u' n _ 1 +f n - 1 u' n _ 1 + u n - 1 u' n _ 1 +r. 


We see that the left-hand side of ( |A8| ) is preserved (up to terms f) if r n _i = (r n — r„_i)/e or 

r n = (r n+ i - r n )/e. 


(AS) 


(A9) 


Note that because of the condition r(0) = 0 the contribution to the Jacobian will be different depending on whether 
we consider positive or negative time. At t < 0 with the regularization rule (|Afl|), Jacobian (A4) equals to 
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J = det 


0 1/e + u' n -1/e 


0 


0 


V.. 


1/e +' 


h n+l 


-l/e 


exp 

/ dtd x u(t, 0) 

/ 



(A10) 


At deriving ( |A1C| ) we have not interested in factors which do not depend on the field u. All the factors should be 
absorbed into the normalization constant Zq introduced in (2.1,2.2). At t > 0 we should go in the positive direction 
in time stating from r(0) = 0. Then we deal with the retarded regularization which does not produce the contribution 
to the Jacobian ( A1C|). With the Jacobian, the action acquires the form 


1 = J dtC, £ = y dx (pd t u + pud x u — upd 2 u) 
+ / dxidx 2 piXi 2 P 2 - i0{-t)d x u(t,0), 


(All) 


where 6 is the step function. 

Now, at calculating the path integral (2.1,2.2) one should use only functions satisfying (|A^) and u(t, 0) = 0. The 
instanton equation for p will acquire the additional term 


d t p + ud x p + vd 2 p — id(—t)5'(x) = S(t)X(x) 


which leads to a nontrivial vacuum: 


ud x p + vd x p 


(A12) 

(A13) 


The requirements that u(x) is a bounded function and p{x) tends to zero w here x —> oo (to provide t he fin ite value 
for c), together with the gauge u(t,0) = 0, determine the unique solution of ( A12 , Alj ). The equation ( A12 ) has the 
first integral: 


2 /2 - v(d x u - a) + c[x(x) - y(0)] = 0 . 


(A14) 


If x —* ±oo then u —> ±rtoo = ±\/2cx(0) — 2ua. The second term under the square root is a small correction of the 
order of Re -1 . For small x , we have u = ax and we find from (A14) the relation a 2 = 2cw 3 . Substituting u —> ax 
into ( A1 3| ) , multiplying it by x and integrating over x, we find 2ac = 1. Then c = l/(2w) and a = ui, that means 
particularly u 00 ~ Llo which is u rms . To obtain the form of p(x), we substitute u —> ax into ( |A1J| ) and solve the 
equation: 


p = I d( exp ' (A15) 

i*i 

The width oi p(x) is L /Re i.e. it is small indeed. The energy of the vacuum is finite E = 3w/4. 

It is important to realize that it is the “gauge” term in the right-hand side of the equation for p that makes vacuum 
non-trivial i.e. not given just by u = 0 and p = 0. In our formalism, the breakdown of Galilean invariance, discussed 
by Polyakov j|], manifests itself due to fixing gauge and making an explicit account of fluctuations of gauge degrees 
of freedom. 

Now, we can repeat all the above instanton consideration with a new vacuum and the additional term in the 
equation for p. That gives no contribution at the main order in large n. 
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